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We present the first determination of charge symmetry violation (CSV) in the spin-dependent 
parton distribution functions of the nucleon. This is done by determining the first two Mellin 
moments of the spin-dependent parton distribution functions of the octet baryons from Nf = 2-1- 
1 lattice simulations. The results are compared with predictions from quark models of nucleon 
structure. We discuss the contribution of partonic spin CSV to the Bjorken sum rule, which is 
important because the CSV contributions represent the only partonic corrections to the Bjorken 
sum rule. 
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I. INTRODUCTION 

Charge symmetry refers to the invariance of the strong 
interaction under a very particular operation in isospin 
space, namely the interchange of u and d quarks and 
also protons and neutrons. Technically, the charge sym- 
metry operator Pes corresponds to a rotation of 180° 
about the 2 axis in isospin space. In nuclear systems, 
charge symmetry is generally valid to substantially bet- 
ter than 1% T]. At the partonic level, charge symmetry 
implies the equality of different parton distribution func- 
tions (PDFs), namely 

uP{x,Q^)^rix,Q^), dPix,Q^)^u"{x,Q^), (1) 

with analogous relations for antiquark PDFs. To date, no 
experimental violation of charge symmetry has been ob- 
served at the partonic level, and current upper limits are 
consistent with the validity of partonic charge symmetry 
in the range 5-10% [21 [3]. 

In this letter we report the first determination of the 
charge symmetry violation (CSV) in the spin-dependent 
parton distributions arising from quark mass differences. 
We begin by extracting the zeroth and first moments 
of the spin-dependent PDFs of the light baryon octet 
by varying the light (degenerate u, d) and strange quark 
masses in a Nf = 2-1-1 lattice simulation. We com- 
pare these results to quark model predictions for the sign 
and magnitude of these moments. Finally, we examine 



the size of the expected contribution of the spin parton 
distributions to the Bjorken Sum Rule. 

II. CHARGE SYMMETRY VIOLATION 

Theoretical models for partonic charge symmetry pre- 
dict that the spin-independent parton CSV distributions 
Su-{x) = uP-{x)~d"'{x) a.ndSd^{x) = dP- {x)~u''^ (x) 
should be roughly equal in magnitude and opposite in 
sign 3J , where the minus superscript denotes the valence 
or C-odd combination of parton distribution functions, 

q^ix) = qix)±qix) . (2) 

The MRST group included valence CSV in a global 
analysis of high energy experimental data [1]. The best 
value obtained in this search was in excellent agreement 
with quark model calculations of valence parton CSV [S] , 
but with very large errors. A recent lattice calculation 
was able to probe the magnitude of CSV violation [5]. 
There, the behaviour of the first moments of the hyperon 
parton distribution functions were studied as the light 
and strange quark masses were varied in a Nf = 2-1-1 
lattice simulation. The first moment of the parton distri- 
butions Su'^{x) and 6d^{x) agreed very well in sign and 
magnitude with both the quark model results and the 
best value from the global fit — with the uncertainties on 
the lattice results substantially smaller than those from 



2 



the global fit. Note that the lattice calculation accessed 
the C-even combinations of partonic CSV distributions, 
so the lattice results contained some sea quark CSV ef- 
fects that were not included in the other investigations. 

We define the m}^ moment of the charge symmetry vi- 
olating spin-dependent quark distributions in the nucleon 
as 



'Au 



/Ad I 



SAd"" = dxx'^iAdPix)- Au'\x)), 



'Ad 



/Am 



(3) 



(4) 



In the limit where the strange and light quarks have 
nearly equal mass, these CSV spin moments are related 
to hyperon spin moments by 



6Au" 
SAd" 



[X- 
(x" 



I An 

\^ 
I As 
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(5) 
(6) 



III. LATTICE SIMULATION DETAILS 



In the numerical calculation of the moments defined in 
Eqs.([3])-([6]), our gauge field configurations have been gen- 
erated with Nf = 2 -f 1 flavours of dynamical fermions, 
using the Symanzik improved gluon action and nonper- 
turbatively C'(a)-improved Wilson fermions [7j. The 
quark masses are chosen by first finding the SU(3)flavour- 
symmetric point where flavour singlet quantities take 
on their physical values and then varying the individ- 
ual quark masses while keeping the singlet quark mass 
Mq = (m„ -I- md + iris) I i = {2mi + ms)/3 constant [Hlini- 
Simulations are performed on lattice volumes of 24^ x 48 
with lattice spacing, a = 0.083(3)fm. A summary of the 
dynamical conflgurations used is given in Table H] More 
details regarding the tuning of the simulation parameters 
can be found in Refs. [SUS]. 

On the lattice we compute moments of the spin- 
dependent quark distribution functions, Aq{x) 



dxx'^iAq^ix) + i-irAq^ix)) , (7) 



where x is the Bjorken scaling variable associated with 
baryon B. This involves calculating the matrix elements 
of local twist-2 operators, namely 

{B{p}\[0'qi'^o...^-^ -Tr]\B{p)) 

-2(x'")ijs{^V^...p^'">-Tr], (8) 



# 






[MeV] 


rriK [MeV] 


1 


0.12083 


0.12104 


460(17) 


401(15) 


2 


0.12090 


0.12090 


423(15) 


423(15) 


3 


0.12095 


0.12080 


395(14) 


438(16) 


4 


0.12100 


0.12070 


360(13) 


451(16) 


5 


0.12104 


0.12062 


334(12) 


463(17) 



TABLE I: Pion and kaon masses on 24 x 48 lattices with 
lattice spacing, a = 0.083(3)fm jSj, where the error on the 
lattice spacing has been included in the errors for and 
rriK- The first column denotes the ensemble number. 



(baryon number conservation). As a result, we only con- 
sidered the first non-trivial moment, {x)q, in our previ- 
ous calculation of the spin- independent CSV |^. The 
lowest moment of the spin-dependent quark distribution 
functions, however, is not protected by such a sum rule. 
Hence, in this work we consider the first two (m = 0, 1) 
moments, which, according to Eq. ([7), contain one C- 
even and one C-odd moment. This a lows us to better 
assess the impact of the sea distribution in our results. 

In this paper we only consider the quark-line connected 
contributions to the first two moments, (1)a(j, {x)Aq, 
which means that we only include the part of coming 
from quark-line connected backward moving quarks, the 
so-called Z-graphs. While the contributions from dis- 
connected insertions are expected to be small (TUl ITT] , 
in the following analysis we will focus on differences of 
baryons and so these contributions will cancel in the 
SU(3)flavour limit and should be negligible for small ex- 
pansions around this limit, as considered here. 

We use a nucleon polarised in the -|-z-direction with 
the standard local operators 



O 



(1) 

Aq 



'^Aq 



and 



O 



(x) 



O 



5{43} 



(9) 



The matrix elements in Eq. 
by considering the ratios: 



' Aq Aq 

are obtained on the lattice 



i?(i,T,p,Og] 



C3pt(t,T,p,oi^>) 



C3pt(i,T,p,0 



i(l)Ag , 



Aq I 



C2ptit,P) 



.mB , . 

~i^{x)Aq- 



(10) 



where C2pt and C^pt are lattice two and three-point func- 
tions, respectively, with total momentum, p, (in our 
simulation we consider only p = 0). The operators 



O^P and O)^' are inserted into the three-point function. 



'Aq 



'Aq 



C3pt{t,T,p,0) at time, r, between the baryon source lo- 
cated at time, t — 0, and sink at time, t. 



IV. LATTICE RESULTS 



where oS^^o-p™ ^ i"^qj^jf^"D - D q . We note that 
in the case of the unpolarised quark distribution func- 
tions the lowest moment is protected by a sum rule 



The operators used for determining the moments of the 
spin-dependent PDFs need to be renormalised, prefer- 
ably using a nonperturbative method such as RI'-MOM 
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FIG. 1: Ratio of doubly-represented quark contributions to a 
baryon's spin, {1)au/{1)au ^^id {1)a3/(1)aii ^ function of 
m^/X^, where the pion masses are normalised to the lattice 
determination of X^. The vertical dotted line indicates the 
physical pion mass. 



# 










1 


1.013(5) 


0.964(22) 


0.958(25) 


1.028(42) 


2 


1.0 


1.0 


1.0 


1.0 


3 


0.999(5) 


1.036(15) 


1.022(30) 


1.009(51) 


4 


0.985(10) 


1.034(26) 


1.054(33) 


0.972(52) 


5 


1.006(16) 


1.064(43) 


1.115(37) 


0.949(66) 



TABLE II: Ratios of the zeroth moment of the C-even, spin- 
dependent hyperon PDFs. 



FIG. 2: Ratio of singly represented quark contributions to a 
baryon's spin, {l)As/{l)Ad W^u/W^d ^ function of 
m^/X^, where the pion masses are normalised to the lattice 
determination of X^. The vertical dotted line indicates the 
physical pion mass. 



# 






{x)as/{x)au 


{x)Au/{x)Ad 


1 


1.036(7) 


0.907(29) 


0.989(28) 


1.025(69) 


2 


1.0 


1.0 


1.0 


1.0 


3 


0.992(7) 


1.064(24) 


1.035(40) 


1.029(79) 


4 


0.956(11) 


1.095(51) 


1.043(34) 


1.000(104) 


5 


0.960(16) 


1.257(93) 


1.067(43) 


0.995(117) 



TABLE III: Ratios of the first moment of the C-odd spin- 
dependent hyperon PDFs. 



[T^l - fH] . Here, however, we will only present results for 
ratios of the first two moments so that the renormali- 
sation constants cancel. In Fig. [T] we present results for 
the ratio of the w(s)-quark contribution to the spin of the 
E(S) baryon to the contribution of the u in the proton, as 



a function of 



normalised with the centre-of-mass of 



the pseudoscalar meson octet, = ^/{2r 



411 MeV. They are also given in Table pl| for each ensem- 
ble. We see that the ratio of the contribution from the 
M-quark to the spin of the E and the proton is roughly 
constant as the quark mass is decreased. We also observe 
that the contribution from the strange quark to the spin 
of the S-baryon is greater than that of the u-quark in the 
proton and increases as the mass of the light (strange) 
quark is decreased (increased). 

Unlike the unpolarised case in [5] , there is no sum rule 
to preserve the total spin-dependent quark contributions. 
This implies that the strange quark contribution to the 
spin of the E doesn't necessarily have to be the same as 
the d-quark to the proton, and in fact we see in Fig. [2] 
that (1)as > (l)Ati- Conversely, we note in Fig. that 
the w-quark in the S baryon feels the effect of the two 
heavier strange quarks and (1)^„ < (l)^^^ decreases as 
we approach the physical point. 



Similar effects are seen in the m = 1 (or x-) moments 
given in Table |III| and shown in Fig. [3] and |4j with the 
exception that in this case, < (a^)Au' 

quark in the E feels the effect of the heavier strange quark 
being present in the baryon. 

To infer the level of CSV relevant to the nucleon, 
we only need to consider a small expansion about the 
SU(3)flavour symmetric point, for which linear flavour ex- 
pansions prove to work extremely well [8]. For instance, 
we can write 



5Au^ = ms 
6Au^ = ms 



dniu 
drriu 



drud 



+ 0{mj), (11) 



where 



i) and we have already made 



use of charge symmetry by equating d{x™')'^^l dm^ 
d{x"')lJdm^ and 9(x'")^^/5m„ = 5(a;")P „/9mrf. 
Similar expressions hold for SA(f and SAd^. 

Near the SU(3)flavour symmetric point, we note that 
the up quark in the proton is equivalent to an up quark 
in a E+ or a strange quark in a 'EP, which we describe 
collectively as the "doubly- represented" quark [15j . 
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FIG. 3: Ratio of doubly represented spin-dependent quark 
momentum fractions, {x)^^/{x)^^ and {x)%s/{xYau ^ ^ 
function of rn^/X^, wiiere the pion masses are normalised 
to the lattice determination of X-j^. The vertical dotted line 
indicates the physical pion mass. 



FIG. 4: Ratio of singly represented spin-dependent quark mo- 
mentum fractions, {x)^^/ {x)^^ and {x)%^/ {xY^^ as a func- 
tion of /X^ , where the pion masses are normalised to the 
lattice determination of X^. The vertical dotted line indicates 
the physical pion mass. 



The local derivatives required for JAu™ can be ob- 
tained by varying the masses of the up and down quarks 
independently. Within the present calculation, we note 



that the difference 



/As 



)au measures precisely 



the variation of the doubly-represented quark matrix ele- 
ment with respect to the doubly-represented quark mass 
(while holding the singly- represented quark mass fixed). 
Similar variations allow us to evaluate the other required 
derivatives, where we write 







'As 


- (x"*) 


p 

Au 


dniu 




rus 


- TO/ 






(x™ 


>Au 


- (X™ 


>Au 


drud 




rris 


- mi 






(x™ 


)au 




>Ad 


drUu 






- mi 








'As 


- (x" 


)P 
'Ad 


drud 




rus 


- mi 





(12) 



With these expressions and Eq. (Ill, wc obtain the rel- 
evant combinations for our determination of CSV in the 
nucleon 



ms- 



)as 



6 Ad"" = ms 



X 



TOs 
m\S + 



TO/ 



TOs - TO/ 



(13) 



By invoking the Gell-Mann-Oakes-Renner relation 
and normalising to the total nucleon isovector m = 0, 1 



spin-dependent moments, we write 
5 Am" ms {{x' 



/An 



/As 



)/ {x"^)Au-Ad 



/ An-Ad 



(to|. - ml)/X^ 



(14) 



SAd" 



ms 



,m\E+ 
As 



i)/ (^™)An-A(i 



/^m\P 
\^ lAu- 



Ad 



[ml mDlXl 



(15) 



Written in this way, the fractional spin-dependent CSV 
terms are just the slopes of the curves shown in Figs. [5] 
and [6] (evaluated at the symmetry point) multiplied by 
the ratio ms/rriq. By fitting the slopes, we obtain 



(5Au" 



ms 



\-^/A«+-Ad+ 


mq 


5Ad°+ 


_ ms 


(-'■)Au+-Arf+ 


mq 


5Av}- 


_ ms 


(^)An--Ad- 


mq 


5Ad^- 


_ ms 


(^)An--Ad- 


mq 



(-0.137 ±0.028) , 



(-0.161 ±0.035) , 



(16) 
(17) 
(18) 
(19) 



Chiral perturbation theory yields the quark mass ratio 
ms/rfiq = 0.066(7) ^16, while the experimentally deter- 
mined moments are (l)An+-Ad+ = 9a = 1.2695(29) [17] 



and (a;)An--Ad- = 0.190(8) [H] in the MS scheme at 
4GeV^. We note that in principle the zeroth moments 
in Eqs. (16) and (TTtI) will receive their scale dependence 



from an addtional term 

z(/i,a) (5Au" + (5AdO) 

3 (l)AM-Aci 



(20) 
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FIG. 5: The zeroth moment of the spin difference between 
doubly and singly represented quarks in the E and H as a 
function of the strange/light quark mass difference. We de- 
duce 5/S.vP and cSAd", respectively, from the slopes of these 
curves (c.f. Eqs. |[l6| - (fl9|). 



where z(ii,a) is the difference between the (scale- 
dependent) singlet and (scale-independent) nonsinglet 
axial-vector current renormalisation constants. At order 
0{a^) in perturbation theory this results in a correction 



of < 1% at 



a 



6 GeV^ [TH [ig and, due to 



its small anomalous dimension, also at other scales, e.g. 
= 4GeV^ 



Substituting these values into Eqs. ( 16 )-( 19 ) yields the 



first lattice QCD estimates of the spin CSV moments 

6Au"+ = -0.0116(27), (5Ad"+ = -0.0036(11) , 
SAu^- = -0.0020(5), 5Ad^- = -0.0009(2) . (21) 

We can make several observations regarding these spin 
CSV moments. First, the fractional spin CSV for both 
moments and both flavours are similar in magnitude and 
all have the same (negative) sign. Second, we can com- 
pare the first moments of the spin CSV distributions with 
the corresponding first moments of the spin-independent 
CSV distributions that were reported in Ref. [6], namely 



6u+ = -0.0023(6), Sd+ = +0.0020(3) , 



(22) 



The first moments of the spin-independent CSV results 
have roughly equal magnitudes but opposite signs, with 
Su being negative and Sd positive, in both qualitative and 
quantitative agreement with quark model predictions [5 , 
[?0j and with the best-fit values from a global fit that 
included valence CSV g]. 

Next we note that the zeroth moments of the spin- 
dependent CSV distributions are larger than the first 
moments. Lastly, we have estimated the CSV associ- 
ated only with the u — d mass difference. It is important 
to also find a method to investigate the CSV induced by 
electromagnetic effects which, at least in the unpolarised 
case, is expected to be of a similar size pTl I22j . 



FIG. 6: The first moment of the spin difference in the E and 
H vs. the strange/light quark mass difference, from which we 
deduce SAu^ and SAd^. 



V. QUARK MODEL COMPUTATION 

We can compare the lattice results with estimates of 
valence quark spin-dependent PDFs obtained from quark 
model calculations. Schreiber, Signal and Thomas [23] 
calculated parton spin distributions from bag models. 
Sather derived an analytic approximation giving va- 
lence parton CSV distributions in terms of derivatives 
of phenomenological PDFs. Sather's equations are valid 
for parton distributions at a low scale appropriate 
for quark models, and should also be valid for CSV spin 
distributions. In this approximation, the valence parton 
CSV spin distributions are 



SAd-{x) = -—-[xAd-{x)] 



SAu'ix) 



M dx 
6M 
Id 



Au^{x) + {l-x) 



dm d 
Af llx 
d 

dx 



Ad-{x), 
Au'{x) 



(23) 



where 6M is the n — p mass difference and Sm is the 
diquark mass difference rridd — rriuu which is determined 
rather accurately to be 4 MeV [2^ . The zeroth moment 
of the spin dependent CSV distributions is overly sen- 
sitive to the small- a; behaviour of these PDFs, a region 
where the quark model results are less reliable. Therefore 
we compare only with the first moments of the CSV spin- 
dependent distributions. Using the model of Schreiber, 
Signal and Thomas we find 

SAu^- = -0.0008, SAd^' = -0.0011. (24) 

Alternatively, if we use the spin-dependent PDFs from a 
Nambu-Jona-Lasinio model calculation [5^1 HZ] together 
with Eqs. (23), we find 



SAu^- 



-0.0003, 5Ad^- 



-0.0007. 



(25) 
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These phenomenological model predictions agree with 
the lattice results insofar as the first spin moments are all 
the same (negative) sign, and have similar magnitudes. 
As for the spin independent case, the result for the down 
distribution, which is dominated by the diquark mass 
shift, is in better agreement than that for the up quark 
where there are a number of small corrections, not all 
included in the Sather approximation. 



possible to observe such a small contribution. However, 
the Bjorken sum rule could in principle be measured at a 
future electron collider, where one could imagine aiming 
for 1% precision [3D]. With such precision it is possible 
that the spin CSV contributions calculated here would 
be sufficiently large to make a measurable difference in 
the sum rule. 

VII. CONCLUSION 



VI. BJORKEN SUM RULE 

Finally, the spin-dependent CSV distributions con- 
tribute to the Bjorken sum rule [28], which has the form 



dx[gl(x) - g-,{x)] 



Ga 



6Gy TT 

Am+ [x) - Ad+ (x) ASAd+ {x) + SAu+ (x) 
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18 



(26) 



In the first line of Eq. ( 26 ) we write the Bjorken sum rule 



in terms of the difference of the spin-dependent structure 
functions 171 for the proton and neutron, integrated over 
all X. In the second line of Eq. ( 26 1 we write the sum rule 



in terms of the first moment of spin-dependent parton 
distributions. This quantity is correct up to terms of 
order 0{as) (there are also higher-twist terms of order 
0(1/(5^)). We see that, except for the CSV corrections, 
this ratio is given by the zeroth moment of the difference 
of the C-even spin distributions Au+ and Arf+ integrated 
over all x. 

We have included the contribution from partonic spin 



CSV in Eq. ( 26 ) which is noteworthy for several reasons. 



First, with the exception of corrections arising from par- 
tonic spin CSV terms, there are essentially no other par- 
tonic corrections to the Bjorken sum rule at leading twist 
(this is one reason why it is so important to obtain precise 
values for this sum rule). Second, the correction involves 
the zeroth moments of 5Am+ and SAd^. At present the 
Bjorken sum rule is best determined from a recent COM- 
PASS experiment at = 3 GeV^ to a precision of about 
8% • Using the zeroth moment obtained from our lat- 
tice calculations (see Eq. (21)) we estimate that the spin 
CSV terms contribute approximately 1% to the Bjorken 
sum rule. At the present measured precision it is not 



In summary, we have performed the first lattice deter- 
minations of the polarised quark moments of the hyper- 
ons, S and S in Nf = 2 + 1 lattice QCD. By examining 
the SU(3)flavour-breaking effects in these momentum frac- 
tions, we are able to extract the first QCD determination 
of the size and sign of charge-symmetry violation in the 
spin-dependent parton distribution functions in the nu- 
cleon, 5Au and SAd. We compare our results with esti- 
mates of the first moment of the parton spin CSV from a 
quark model calculation, obtaining qualitative agreement 
with the quark model results. Finally, we estimate the 
contribution of partonic spin CSV to the Bjorken sum 
rule, and show that spin CSV effects should change the 
Bjorken sum rule by approximately 1%. 
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